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Abstract 

The branching ratios of the r~ — ► h~ u T , t~ — > h~ ir°v T and r~ — > h~ > 2tt°u t decays have been measured 
using the 1991-1995 data recorded with the OPAL detector at LEP. These branching ratios are measured 
simultaneously using three selection criteria and are found to be 

BR(t-^Ti/ t ) = (11.98 ±0.13 ±0.16)% 

BR(r~ — > h~ vr°z^ T ) = (25.89 ±0.17 ±0.29)% 



BR(t~ -> h~ > 2vr V) = (9.91 ± 0.31 ± 0.27) % 



o 



where the first error is statistical and the second is systematic. 
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1 Introduction 



Precise tests of the Standard Model can be made using r leptons Measurements of the r allow the 
study of the structure of the weak currents and the universality of the couplings of the charged leptons 
to the gauge bosons, and can be used to search for evidence of new physics. Further, the r is the only 
lepton to decay into hadrons, allowing the study of the strong interaction. For many of these studies an 
accurate knowledge of the properties of the r is essential. In this paper we report on a measurement of the 
branching ratios of the r~ —* v T , t~ — > h~ ir°v T and r~ — » h~ > 2ir°v T decays [] using the 1991-1995 
data recorded with the OPAL detector at LEP. 

The measurement of the one-prong hadronic branching ratios is made by selecting a sample of tau decays 
with one track (one-prong) and then counting the number of 7r°'s in each decay. The one-prong decays 
were subdivided into samples with 0, 1 or > 2 7r°'s from which the branching ratios for the three signal 
channels were determined. In Table |l| we list the r decays that are included in each signal channel. We 
use the PDG definitions for these signal channels [0]. Note that not all modes were included in the Monte 
Carlo simulation. No attempt to separate charged pions and kaons was made in this measurement. The 
limited granularity of the OPAL electromagnetic calorimeter made it impossible to resolve unambiguously 
the decays with two 7r°'s from those with three or more 7r°'s with the 7r° identification algorithm used in 
this work, hence only a measurement of the > 2tt° branching ratio is presented. 



Selection 




Decay 


mode 


Weight 


Comment 


t- -► hr 


V T 


T~ — ► 












T~ — ► 


K _ i/ T 






r- -► hr 




r~ — > 


7T _ 7T U I/ T 










r~ — ► 








t- -» hr 


> 2vrV 


r~ — * 


7T~ 27T U V t 










T~ — ► 


K~ 2vr° z/ T 










r~ — > 


7T~ 37T° f T 










r~ — > 


K~ 3vr° zv r 




Not modelled 






r _ — > 


ft - 4vr° 




Not modelled 






r~ — > 


-x7° 


0.157+ 








r~ — > 


k~kV 


0.157+ 


Not modelled 






r~ — ► 


-t7° 


0.157+ 








r~ — > 


K-K°7r°^ r 


0.157+ 


Not modelled 






r~ — > 


vr-K K°^ r 


0.0246+ 








r~ — > 


7T~7T ?7 Z/r 


0.319* 





f Only the Kg — > 2ix decay included. 
* Only the r\ — * Sir decay included. 



Table 1: The various decay modes for each selection. The weights reflect the fraction of each decay that 
contributes to the signal. The weight is equal to unity when no number is explicitly given. A number of 
the decay modes were not simulated in the Monte Carlo and these are indicated in the last column. 



Charge conjugation is implied throughout this paper. The symbol h is used to indicate either it or K . 
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2 OPAL detector 



A detailed description of the OPAL detector can be found in Ref. ||. A description of the features relevant 
for this analysis follows. 

A high precision silicon microvertex detector surrounds the beam pipe. This covers an angular region of 
| cos#| < 0.8 and provides hit information in the x — y (and z after 1992) directions^] Q. Charged particles 
are tracked in a central detector enclosed inside a solenoid that provides a uniform axial magnetic field of 
0.435T. The central detector consists of three drift chambers: a high resolution vertex detector, a large 
volume drift jet chamber and the z-chambers. The jet chamber records the momentum and energy loss of 
charged particles over 98% of the solid angle and the z-chambers are used to improve the track position 
measurement in the z direction ||. 

Outside the solenoid coil are scintillation counters which measure the time-of-flight from the interaction 
region and aid in the rejection of cosmic events. Next is the electromagnetic calorimeter (ECAL) that is 
divided into barrel (|cos#| < 0.82) and end-cap (0.81 < |cos#| < 0.98) sections. The barrel section is 
composed of 9440 lead-glass blocks each subtending approximately 40 x 40mrad 2 and with a depth of 24.6 
radiation lengths. 

Beyond the electromagnetic calorimeter the iron of the solenoid return yoke is segmented into layers and 
instrumented with limited streamer tubes as the hadron calorimeter (HCAL). In the region | cos#| < 0.81 
this detector typically has a depth of 8 interaction lengths. Beyond the hadron calorimeter is the muon 
chamber system, composed of four layers of drift chambers in the barrel region. 



3 Event selection 

The results presented in this paper are based on the data taken during the 1991-95 runs with the OPAL 
detector at LEP. Approximately 90% of the data were taken at a centre-of-mass energy equal to the mass 
of the Z° boson Mz, with the remaining data taken within 3 GeV of Mz- 

The Monte Carlo samples used in this analysis consist of 300 000 r-pair events generated at y/s = Mz 
and two samples of 100 000 events each generated respectively at 2 GeV above and below Mz- The Monte 
Carlo samples were generated with KORALZ 4.0 || and TAUOLA 2.0 J5J and then processed through the 
GEANT H OPAL detector simulation Q. For this analysis an admixture of events from the Monte Carlo 
samples generated above and below Mz are added to the events generated at Mz to reflect the distribution 
of centre-of-mass energies in the 1991-95 data set. 

It has been found in a recent analysis []l0| that the model of Kiihn and Santamaria (KS) pT| , which 
is used in TAUOLA 2.0, does not satisfactorily describe the dynamics of the tau decay through the a\ 
resonance. The model of Isgur et al. (IMR) fl2|| , although also not providing a completely satisfactory 
description of the r~ — > a±v T decay, was found in [|hJ to be in better agreement with the data than the KS 
model. This improvement was found to be due mainly to the inclusion in the IMR model of a polynomial 
background term, which accounts for (13.8 ± 2.4)% of the total decay. 

For this analysis, we have therefore chosen to apply weights to the Monte Carlo generated events so 
that they are distributed dynamically like the IMR model description. The weights applied are based on 

2 The OPAL coordinate system defines the +z axis in the e~ beam direction. The angle 8 is measured from the +z axis 
and <f> is measured about the z axis from the +x axis which points to the centre of the LEP ring. 
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the values of Q 2 , s\, and S2 at the tree level, where Q 2 is the invariant mass squared of the 7r~ 2n° system. 
The Dalitz plot variables si and S2 are defined in terms of the pion 4-momenta as si = (p2 + P3) 2 and 
s 2 = (Pi +Pz) 2 -, with the labels chosen such that p% refers to the charged pion. The weighting has been 
accomplished by generating normalized three dimensional arrays of decay rate T(Q 2 , s\, S2) for each of the 
two models, and taking the bin-by-bin quotient as the weight factor. The array grid size was chosen to 
be 0.05 GeV 2 for each of Q 2 , s±, and S2- The model parameters were taken from [10] for the IMR model, 



while KS model parameters, taken from were those used in TAUOLA 2.0. A smoothing algorithm 
was applied to the IMR array to force the polynomial background term to zero smoothly at the s± and 
S2 physical boundaries, but in such a way as to preserve the measured Q 2 dependent shape and fractional 
contribution of the polynomial background term [|lC|]. 

The event selection starts by identifying e + e~ — > t + t~ events, then r decays with one track are 
identified and clusters are formed in the electromagnetic calorimeter. The 7r 0, s in each one-prong jet are 
reconstructed and background r decays are rejected. 



3.1 Tau pair selection 

The r-pair sample is created by selecting events in the barrel region of the detector with two back-to-back 
cones or jets [ 13 1 , of half-angle 35° degrees. Each event is required to have 



-EECAL + -Etrack > 0.01 EqM 



| cos 6»| < 0.68, 

where -Btrack is the scalar sum of the momenta of the tracks, £ecal is the total energy of the clusters in 
the electromagnetic calorimeter, Eqm is the centre-of-mass energy of the e + e~ beams and |cos#| is the 
average value of | cos 9\ for the two jets. 

Cosmic and beam related backgrounds are rejected by placing requirements on the time-of-flight detec- 
tor. Additional requirements are needed to separate the t + t~ events from other two fermion background 
(e + e~ — > ff) events: 



• Multihadronic events (e + e~ — ► qq) at the LEP energies are characterized by large track and cluster 
multiplicities. These events are rejected by requiring at least two and not more than six tracks and 
not more than ten electromagnetic calorimeter clusters. 

• Bhabha events (e + e~ — * e + e~) are characterized by two back-to-back high energy charged particles 
that deposit close to centre-of-mass energy in the electromagnetic calorimeter. Bhabha events are 
rejected by requiring the r-pair candidates to have -Eecal < 0.8-Ecm or -Eecal + 0.3-E t rack < Eqm- 

• Muon pair events (e + e~ — ► are identified as two high momentum back-to-back tracks that 
leave little energy in the electromagnetic calorimeter. These events are removed if the tracks have 
associated activity in the muon detectors or hadronic calorimeter and £ecal + -^track > 0.6Eqm- 



In addition to the background from two fermion events, there are also two-photon events (e + e~ — > 
(e + e~)A, where X = e + e~, t + t~, qq) that must be rejected. Two photon events leave little 

energy in the detector as the e + and e~ particles are emitted at angles close to the beam and are often 
undetected. In addition, the detected particles tend to have a large acollinearity angle. The acollinearity 
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Background 


Contamination (%) 


e + e~ - 


-» ^ n 


0.72 ±0.05 


e + e~ - 


-> e + e~ 


0.41 ± 0.07 


e + e~ - 


-> qq 


0.28 ± 0.04 


e + e~ - 


-> (e + e~)e + e~ 


0.07 ±0.02 


e + e~ - 


* (e + e~)/i"V~ 


0.08 ± 0.02 


Total 




1.56 ±0.10 



Table 2: Non-tau background in the r-pair sample. 



angle, 6 aco i, is denned to be the complement of the angle between the two tau jets in the event. These 
events are rejected by requiring 

#acol < 15° , 
E vis > 0.03-Ecm j 

where E v i s is the sum of the visible energies of the jets (taken for each jet as the maximum of the sum of 
the track and electromagnetic calorimeter cluster energies). If E v [ s < 0.20i?cM then events are rejected if 
they satisfy 

Pj acks > 2.0GeVor Pj CAL > 2.0 GeV , 

where -bracks (-Pecal) i s the vec tor sum °f momentum (energy) of all tracks (electromagnetic calorimeter 
clusters) in the transverse direction. 

The tau selection applied to all data yields 95364 r-pair candidates. The non-tau background contri- 



butions in the sample have been investigated in reference [14] and are shown in Table The e + e~ — > qq 
background has been re-evaluated by comparing the number of clusters in the electromagnetic calorimeter 
with those predicted by the r Monte Carlo. The e + e~ — > qq events that enter the r sample tend to have 
a larger number of clusters than r-pair events. 



3.2 One-prong selection 



Jets with one track are identified as one-prong decays. Jets with 2 or 3 tracks may also be identified as 
one-prong decays if the additional tracks are associated with a photon conversion. The track not identified 
as a conversion electron is henceforth called the primary track. The photon conversion algorithm used in 
this analysis is described in [15]. A total of 158971 one prong jets are selected of which 5960 are jets with 
2 tracks and 1903 are jets with 3 tracks. 



3.3 Clustering algorithm 



A fine clustering algorithm [16] is used to identify the particles in the one-prong r decays. The fine 
clustering algorithm limits the cluster size in the electromagnetic calorimeter to be 2 x 2 blocks in 9 and 
<p. Clusters adjacent to other clusters may have fewer blocks. Both data and Monte Carlo show that, on 
average, 99% of the energy of an electron and 95% of the energy of a charged pion that is deposited in the 
lead-glass calorimeter is contained in the 2x2 cluster. 

The electromagnetic calorimeter clusters are matched to the tracks using a x 2 significance parameter in 
6 and <f> that is weighted by the uncertainty in the track position and the uncertainty in the cluster centroid. 
An electromagnetic calorimeter cluster that is not associated to a track and that has energy (Snc) greater 
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than 0.5 GeV is classified as a neutral cluster. The energy of the Monte Carlo electromagnetic calorimeter 
clusters has been smeared so that the energy resolution in the data and Monte Carlo are approximately 
equal. Figure [l] shows the distribution of the number of neutral clusters per one-prong jet ^. 

In about 1% of r decays there is a neutral cluster in the jet that is created by a radiative photon. If 
the invariant mass of the primary track and any cluster is greater than 3.0 GeV, then that cluster is not 
considered a neutral cluster. Figure [2] shows the distribution of the energy of neutral clusters removed by 
this requirement. 

3.4 7T° identification 

The 7T° decays 98.8% of the time into two photons with the remaining 1.2% of the decays into an e + e~7 
final state (Dalitz decay). The algorithm used to identify 7r°'s is applied in four sequential steps: 

1. Any neutral cluster in the jet with E^c > 9.0 GeV is identified as a 7r°. In the selected samples of 
t~ — > h~ tt°u t and r~ — > h~ > 2tt°u t jets, approximately 55% and 34% of the 7r°'s, respectively, are 
identified by this criterion. 

2. Pairs of neutral clusters, or a neutral cluster and photon conversion, each with neutral cluster or 
photon conversion energy less than 9.0 GeV, are candidates to form a ir°. The pair is considered a 
7T° candidate if its energy is at least 3 GeV and its invariant mass (m 77 ) is consistent with the 7r° 
mass using a x^o requirement. The x\a variable is defined as 

2 = (m 77 - m n o) 2 

Xti-0 — 2 ' 

where m n o is the mass of the tt° meson and cr m77 is the calculated uncertainty on m 77 . A pair is 
considered a tt° if x^o < 9.0. 

The number of 7r°'s formed in each jet using this method is not limited. If there is an ambiguity 
between neutral clusters and photon conversion, the combination that gives the best x^o is chosen. 
Figure || shows the invariant mass distribution of two neutral clusters before and after this selec- 
tion criterion is applied. In the selected samples of r~ — ► h~ and t — > h > 2ir®v T jets, 
approximately 17% and 14% of the 7r°'s, respectively, are identified by this criterion. 

3. Any remaining neutral clusters with -Enc > 2.2 GeV are classified as 7r°'s. In the selected samples 
of — * h~ and t — ► h > 27r°v T jets, approximately 27% and 33% of the 7r 0, s, respectively, 
are identified by this criterion. 

4. Frequently the 7r° cannot be resolved from a track. If the cluster associated to the track satisfies 
both 

E/p > 1.0, 

and E - 0.3p > 2.2 GeV , 

then we consider the cluster to be an overlap 7r°. The energy of the 7T° iS estimated to be E — 0.3p 
where E is the energy of the cluster (charged hadron plus 7r°). The energy deposited by the charged 
hadron in the electromagnetic calorimeter is on average one-third of its momentum (0.3p). These 

The Monte Carlo distributions shown in all figures are normalized to the number of events after the r-pair selection. The 
t branching ratios used are those published by the PDG Q except for the signal channels where the results of this analysis 
are used. 
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7r°'s are only permitted in jets where one or more 7r°'s are identified by any of cases (1) to (3) above. 
If one or more conversion tracks point to the electromagnetic calorimeter cluster associated to the 
primary track then the cluster energy (E) is modified by subtracting the momenta of the conversion 
track(s). In the selected samples of r~ — » h~ > 2-k°v t jets, approximately 19% of the 7r°'s are 
identified by this criterion. 

Each one-prong jet is classified as either r~ — > h~~ v T , t~ —> h~ tt q v t or - -> h~ > 2ttV- Fi gure ffl 
shows the distributions of the number of 7r°'s per jet and the types of 7r°'s identified by this algorithm 
after background rejection (described below). The energy distribution of 7r°'s in r~ and 
T~ —> h~ > 2ir°v T jets is shown in figure || 



3.5 Background rejection 



A number of additional requirements are applied to remove residual backgrounds. In particular, the 
r~ — > h~ u T sample has contamination from r~ — * e~ V e v T and r~ — * \i~ v T decays. This contamination 
is reduced by requiring 

E/p < 0.75, 
P/E heam > 0.05, 

A rlayers _ n 
JV MB _ u ' 

where -Ebeam is the energy of the LEP beam and N^^ rh is the number of layers hit in the muon barrel 
detector. In addition, backgrounds from e + e~ — ► e + e~ and e + e~ — ► events in the t~ — > h~ v T 

sample are rejected by removing events where the acoplanarity angle between the two jets in the event is 
less than 0.003 radians and the primary tracks in each jet have p > 30 GeV. The acoplanarity angle is 
defined to be the complement of the angle between the two jets in the transverse plane of the event. 

In the t~ — > h~ and t — > h > 2ir°v T samples the main backgrounds are due to misidentification 
of other signal channels. In the t~ — > h~ and the r — ► h > 2tt u t samples, the background is 

reduced by requiring the invariant jet mass to be less than 2.0 GeV. In the t~ — > h~ Tt°v T sample, the 
invariant jet mass is also required to be greater than 0.4 GeV. The invariant mass distributions of the 
r~ — > hr 7T°^ T and —> h~ > 2is®v T samples are shown in figure ^. 

A number of r _ — ► K decays are mistakenly selected into the r — > h > 27r°u T sample. The 

total energy of the 7r°'s from the r~ — > h~ ifiv T decays is relatively low and this contamination is reduced 
by requiring 

prtotal 

> 0.8, 

p 

where -E*° tal is the energy sum of all 7r°'s identified in the jet and p is the momentum of the primary track. 



4 Estimation of r backgrounds 



The backgrounds from t~ — > e~ 17 c u T and r~ — > pT v T decays in each sample are measured using 
distributions previously unused in the selection process. A region in each distribution that is dominated 
by the background is used to make the measurement. The data and Monte Carlo distributions in the 
background-dominated region are compared and any deviations between the two are assumed to be caused 
by the background in question. If the ratio of data to Monte Carlo is different from unity then the Monte 
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Selection 


Background 






r — > e f e v r 


T~ -> M~ ^ 




— > h~ v T 


1.30 ±0.18 


1.39 ±0.26 


T~ 


— ► h~ TT°U T 


1.04 ±0.06 


0.91 ±0.28 


T~ 


-> h~ > 2tt°u t 


0.80 ±0.28 





Table 3: Correction factors used to scale the r — > e z^ e ^ T and r — > v T backgrounds. 



Background 




Selections 






r ^ h v T 


r~ — > h~ 7r u u T t~ 


-> h~ > 2vrV 


t~ — * e~ V e v T 


1.43 ± 0.20% 


1.59 ±0.11% 


1.67 ±0.59% 


T~ -> //~ zv r 


1.57 ± 0.29% 


0.68 ±0.21% 


0.02 ± 0.04% 


T - -> /T K° X zy r 


1.80 ±0.14% 


1.57 ±0.11% 


1.81 ±0.15% 


Other background 


0.07 ± 0.02% 


0.53 ± 0.03% 


3.32 ± 0.24% 


Total 


4.88 ± 0.38% 


4.37 ± 0.26% 


6.82 ± 0.65% 



Table 4: The corrected backgrounds as a percentage of each selection. 



Carlo predicted r background is rescaled in the background calculation. The error of this correction factor 
is the combined statistical errors of the data and Monte Carlo. The uncertainties on the correction factors 
are included in the background errors (including the case where the correction factor is unity). 

The r~ — ► e~V c v T background is estimated using the energy loss (dE/dx) in the jet drift chamber. 
Figure [7| shows the dE/dx distribution in the data and Monte Carlo for jets with one ir° and p < 5 GeV. 
The ratios of data to Monte Carlo r~ — ► e~ V c v T background jets are calculated using jets that have 
9.25 < dE/dx < 12.0 keV/cm. The ratio obtained in each sample is used as a correction factor and is 
shown in Table ^. 

The t~ — > pT Vn v T background in the r~ — ► h~ u T sample is measured by creating a sample of 
r~ —* pT Vn, u T jets in both data and Monte Carlo. The efficiency of the muon chamber requirements 
to reject r~ — * p~ v T is measured using these samples and the ratio of the efficiencies is used as the 
Monte Carlo scaling factor for r~ — > p~ V^ v T background in the r~ — > h~ v T sample. The r~ — > p~ Vn v T 
contamination in the r~ — ► h~ sample is measured by identifying tracks with hits in the muon 

chambers, since the muon chambers were not used in this selection. The ratio of data to Monte Carlo is 
used as the scaling factor for r~ — > pT V^ v T background in the r~ — * h~ sample. 

Table ^| shows the estimated backgrounds in each selection. In addition to the t~ — > e~ V e v T and 
r _ — ► p~~ Vfj, v T decays, there is background from r~ — > h~ K X v T decays where X is any number of neutral 
particles in the final state. There is also background from other tau decays, labelled as "Other background" 
in Table [|, which are mainly tau decays to r\ (27%) and u modes (60%), with lesser contributions coming 
from decays with three charged hadrons in the final state (13%). Both the r~ — > h~ K X&v and other tau 
decays are measured from Monte Carlo information. The errors are calculated from the uncertainties in 
the efficiency matrix and the background branching ratio from the PDG [||. The background uncertainties 
in each selection, shown in Table ^, are used in calculating the systematic errors. An additional modelling 
uncertainty, described in Section]^, is added to the r~ — > /i~K°Xi/ T decays. 
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Selection Selection efficiency from MC 

t~ — > h~ v r t~ — > h~ t^v t t~ — > h~ > 2tt^u t 
Ovr u 0.5502 ± 0.0024 0.1003 ± 0.0010 0.0249 ± 0.0008 

lvr 0.0455 ± 0.0010 0.5985 ± 0.0016 0.4273 ± 0.0025 

> 2tt 0.0005 ± 0.0001 0.0237 ± 0.0005 0.2730 ± 0.0023 

Table 5: Efficiencies for identifying the signals for each selection. The errors on these efficiencies are based 
on Monte Carlo statistics only. 



5 Branching ratio calculation 



The branching ratios are calculated using the information from each selection simultaneously. Each selec- 
tion can be expressed in terms of the efficiency for detecting each decay mode, the branching ratio of each 
mode and the number of events selected in the data. For each selection i, the equation is written as 

M Ni 
eu-Bi + e i2 B 2 + e i3 B 3 + J2 e ik B k = _ j non _ T ^ (1) 

where (j = 1,3) are the efficiencies for selecting signal j using selection i and (k = 4, ...) are the 
efficiencies for selecting the r background modes using selection i. Bj (j = 1,3) are the branching ratios 
of the signal channels and B^ {k = 4, ...) are the branching ratios of the backgrounds. Ni is the number 
of data events that pass the selection i, j non ~ T is the fraction of non-tau events in the r-pair sample and 
N T is the total number of data r's that pass the r-pair selection. The selection efficiencies (ey) for both 
signal and background are determined from Monte Carlo with the main backgrounds corrected using data 
distributions (described in section |3|). The r background branching ratios are taken from the PDG and 

pon-r is ghown in Table U 

The solution of the three simultaneous equations gives the branching ratios in the tau-selected sample. 
These branching ratios are then corrected to account for the biases introduced to the r-pair sample by 
the r-pair selection. These factors are 0.989 ± 0.002, 1.019 ± 0.001 and 1.013 ± 0.002 for the r~ -> h~ v T 
t~ — > h~ ir v T and r _ — > h~ > 27r°v T decays, respectively. 



6 Results 



The r-pair selection identifies 190728 r candidates. The selections for r~ — > u T , t~ — > h~ tx^v t and 
r _ —> h~ > 2ir®v T yield 18547, 40537 and 6802 events respectively. The efficiencies for detecting these 
signals are listed in Table [|. The efficiency for detecting r~ — > h~ v T is an average of the efficiencies for 
selecting r~ — ► ix~ v T and r~ — > K~ v T weighted by their relative branching ratios. Similarly, the efficiency 
for detecting r~ — ► h~ ir°v T is an average of the efficiencies for selecting r _ — ► ir~ ir°v T and r~ — > K" 
weighted by their relative branching ratios. The same method is used for the r~ — ► h~ > 2ir°v T efficiency 
where the decay modes are listed in Table [l|. 

The backgrounds in the r _ — ► h~ u T , t~ — > h~ ir v T and r~ — > h~ > 2it v t samples are given in 
Table ||. The branching ratios for r~ — > h~ u T , r~ — > h~ n®v T , and r~ — > h _ > 2-7r i/ r are calculated 
using the method described in section || and the results are shown in Table ||. The measurements of these 
branching ratios are correlated. The correlation coefficients between branching ratios calculated using the 
statistical errors on each branching ratio are given in Table ^. 
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Channel 




Branching Ratio (%) 


r~ -» /i" 




11.98 ±0.13 ±0.16 


r - -» /r 




25.89 ±0.17 ±0.29 




> 2ttV 


9.91 ±0.31 ±0.27 



Table 6: Branching ratio results. The first error is statistical and the second systematic. 



Sample 


r — ► /i 7r u z/ T r 


/T" > 2tt u i/ t 


r~ — * h~ u T 


-0.430 


0.167 






-0.470 



Table 7: The correlation coefficients between each measurement. 



7 Systematic errors 

The systematic errors on the branching ratios are shown in Table ||[ The efficiencies for detecting the 
signals and backgrounds are taken from Monte Carlo and special attention is paid to the Monte Carlo 
modelling. The first half of Table || gives the systematic errors due to aspects of the analysis such as Monte 
Carlo modelling of hadronic showers, track finding, conversion finding, radiative photons and energy scale. 
The second half of Table || gives the systematic errors due to the backgrounds. We briefly describe the 
individual contributions to the systematic errors. 

The a\ — > 7r~2-7r decay was simulated using the IMR model, as described in section ||. The IMR model 
gives a better description of the — > it~it + it~ data primarily because of the inclusion of the polynomial 



background term [12]. The uncertainty in the contribution of the polynomial background was measured 
in |l(| to be approximately 17%. The systematic error in the branching ratios due to the uncertainty in 
the a% model is conservatively estimated by varying the normalization of the polynomial background by 
±50% (3a). 

The Monte Carlo models the hadronic showering in r decays reasonably well. The thresholds for 
the neutral clusters, one-cluster and two-cluster 7r°'s, as well as other requirements, have been carefully 
studied to avoid any poorly modelled regions due to low energy hadronic clusters. However, as the energy 
thresholds are lowered, deviations between the data and Monte Carlo appear since clusters close to the 
track from hadronic interactions are accepted into the 7r°-finding algorithm. The lower energy thresholds 
of both the one and two cluster 7r° cases are varied to estimate the systematic error due to the Monte 
Carlo modelling of energy deposition in the electromagnetic calorimeter. The lower energy threshold of 
the one-cluster 7r° case is varied from 2.0 to 3.0 GeV and the energy threshold of the two-cluster 7r° case 
is varied from 2.6 to 4.0 GeV. The maximum variation from each nominal branching ratio is taken as the 
systematic error in each case. 

The effect of tracks passing close to the anode plane of the OPAL jet chamber is considered as a source 
of systematic error. The Monte Carlo does not perfectly model the position of tracks close to the anode 
plane. The branching ratios were recalculated by removing tracks within 0.25° of the anode plane. The 
resulting change in the branching ratios is taken as the systematic error. 

Approximately 5% of the identified 7r°'s are composed of a neutral cluster and a photon conversion. 
Although the Monte Carlo does a reasonable job of modelling conversions, there are minor discrepancies 
between data and Monte Carlo. The ratio of data to Monte Carlo jets with an identified conversion pair 
is found to be 0.985 ± 0.013. The branching ratios are calculated with Monte Carlo jets containing a 
conversion weighted to simulate a ±3<t (±4%) change in conversion identification efficiency. The maximum 
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Systematic error for each selection (%) 




T — > h ~U T 


T — > h TT^V 

1 ' 1 ll 1 \ IS rj- 


r~ — > h~ > 2n IJ v 

1 'IV ■ — £u / i 


MC statistics 


n 1 1 


1 7 


14 


n i m on pi 1 1 n c 


n m 

u.ux 


1 3 




1-rliiRtpr 7r^ t,Ti TP-fili nl H 


0.03 


0.02 


0.03 


2-cluster 7r^ threshold. 


0.03 


0.11 


0.11 


Anode plane cut 


o on 




03 


T-* not on poiTVPvmoTm 


0.06 


0.03 


0.02 


Energy smearing 


0.01 


0.01 


0.03 


Radiative clusters 


0.03 


0.00 


0.00 


Energy scale 


0.03 


0.03 


0.04 


Signal BR uncertainty 


0.01 


0.05 


0.03 


Bias factors 


0.02 


0.03 


0.02 


Unmodelled channels 


0.00 


0.02 


0.04 


Non-r backgrounds 


0.02 


0.04 


0.01 


r~ — > e~ Z? e 


0.04 


0.07 


0.08 


T~ ->■ M~ F M 


0.05 


0.08 


0.01 


r - -» K° X v T 


0.06 


0.07 


0.04 


Other backgrounds 


0.00 


0.03 


0.03 


Total 


0.16 


0.29 


0.27 



Table 8: Systematic errors. 



variation from each nominal branching ratio is taken as the systematic error. The sensitivity of the results 
to photon conversions is also checked by dropping the conversion routine and recalculating the branching 
ratios. 

The energy of each electromagnetic calorimeter cluster in the Monte Carlo is smeared. The uncertainty 
due to this smearing is assessed by varying the amount of smearing applied by ±20%. The change in the 
branching ratios is taken as the systematic error. 

Clusters were considered to be due to radiative photons if the mass of the track and cluster was greater 
than 3 GeV. These clusters are ignored by the 7r°-algorithm. This requirement was removed and the change 
in the branching ratios was taken as the systematic error. 

The energy scale uncertainty reflects the uncertainty in the electromagnetic calorimeter calibration 
between Monte Carlo and data. To determine the systematic uncertainty the clusters are rescaled by 
±0.2% and the largest effect on the branching ratios is taken as the systematic uncertainty. 

The selection efficiencies depend on the relative branching ratios of the individual decay modes. The 
signal branching ratios were calculated using the PDG branching ratios and uncertainties 0. The depen- 
dence of the efficiency and the decay branching ratios on the signal branching ratios was evaluated and 
included as a systematic uncertainty. 

A small correction must be applied to the branching ratios to correct the slight bias introduced by the 
r-pair selection criteria. The dependence of the bias factor on the r-pair selection and on the branching 
ratios of each channel is found to be relatively small. In addition, other Monte Carlo samples using a 
different electromagnetic shower model give similar results. The systematic error on each branching ratio 
is calculated directly using the bias factor error. 

The Monte Carlo used in this analysis did not include some decay modes defined as signals by the 
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PDG |^| . Table |T| shows these decay modes. The systematic effect of the t~ — > h~ ATr°v T mode on the 
r~ — > /i~ 7r°^ r and r _ — > /i~ > 27r°z/ T samples is calculated using the PDG branching ratio and the 
assumption that the efficiency for r~ — ► h~ £ir°i/ T is the same as that for r~ — > 3tt°v t . The uncertainty 
for each of these modes is assumed to be 50% of their respective branching ratios. The total systematic 
error is calculated from these uncertainties. 

The systematic errors on the branching ratios due to non-r and r backgrounds are calculated from 
the uncertainties given in Tables || and [|, respectively. An additional error has been added to the r~ — > 
h~ K° X v T decays to account for any uncertainty in the energy deposited by a neutral hadron in the 
electromagnetic calorimeter. A 10% error was added to the r~ — > h~ K Xz/ T background in the t~ — > z/ T 
sample. A second error was added to the t~ — ► h~ K° X z/,- background in all three samples to account for 
possible migration of the background from one ir° sample to another. 



8 Discussion 



The r — > h v T branching ratio is measured to be (11.98 ± 0.13 ± 0.16) % and is compared with other 
published results in figure ||. The r~ — > h~ v T branching ratio agrees well with the PDG world average of 
(11.70 ± 0.11)% J§ but is slightly above the recent CLEO measurement of (11.52 ± 0.13)% g7) and the 
theoretical prediction of (11.65 ± 0.06)% made by Decker and Finkemeier (H 



The r~ — > h~ v T branching ratio can be used to measure the ratio of the charged current coupling 
constants of taus and muons. Lepton universality requires that the weak charged current gauge coupling 
strengths be identical: g e = = g T . The g^/g T ratio is probed by comparing the decays h~ — ► li'v^ and 
r _ — > h~ v T . An expression for g^/g^ is given by 



2ml \ BR(r _ ^h~v T ) 



mi 



where 



(l + S h 



f T T m T \ 
\Thm h J 



1 - {m h /m r y 



BR (h -»• ^ M ) , 



"^u and K 



1 - (m^/rrih) 

and 5ft is an electromagnetic radiative correction |^8[ . The tt~ — > [i 

the masses of the tau, muon and the pion, and the lifetime of the pion are taken from the PDG B 
ratio g T /g^i is found to be 

1.018 ±0.010, 



„ branching ratios, 
The 



9t_ 

9im 



where the error on g T / 'g„ is dominated by the uncertainties in the r~ — > /i~ z/ r branching ratio and the 
OPAL t lifetime p9| . The point in figure |9| shows the result from this work plotted against the r lifetime. 
The Standard Model prediction is shown as the shaded band with a width reflecting the uncertainty in 
m T . If we recalculate the world average t~ — ► h~ v T branching ratio including our result and use the PDG 
value for the tau lifetime @, then the ratio g T /g^ is found to be 1.003 ± 0.005. 



v c Va decays. 



The ratio of g T /g^ can also be measured by comparing the r~ — * e~V c v T and fi — * 
Using the PDG values for the r~ — > e~T> e v T and r lifetime gives g T /g^ = 1-000 ± 0.003. Although this 
result is more precise, the two measurements are complementary. The measurement using the r~ — * n~ v T 



and 7r 



H v. 



decays probes the couplings to a longitudinal W boson, while the measurement using the 



v ZA 



and /j, 



e i/ e ^ 



decays probes the couplings to a transverse W. 
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The decay t~ — > /i~i/ T is also sensitive to the presence of physics beyond the minimal Standard Model. 
As an example, in i?-parity violating (RPV) extensions of supersymmetric models [2C] r~ — > h~v T and 
h~ — > n~Vn receive an additional contribution from the exchange of a right-handed scalar d-type quark, 

A form leading to a modification 
^) as outlined in plfl . 



d^, where A; is a family index. This exchange can be cast in a V — A® V 
of the Standard Model predicted couplings by a term proportional to 

The RPV Yukawa couplings probed by t~ — ► h~v T and h~ — > decays are A 31fc and A 21fc respectively. 

Note that, under the assumption that only one RPV Yukawa coupling is non-zero, the decays r~ — ► K~z/ r 
and K~ — ► yTVa are unaffected by RPV since they involve products of two different coupling constants. 



Using the formalism developed in El], limits can be set on A 31fc and A 21fc from the expression, derived 
here for the first time, 



\\r 12 _ I \/ \2 

' 31fc ' ' f fc ' = 2v / 2G F y ud 
m A (d K R ) 



2m M 2 BR(r" -> /jTi/ t ) 



1 



Making the assumption that only one Yukawa coupling is non-zero leads to the following 95% confidence 
level limits, calculated using m(d R ) = 100 GeV, of 



A' 31fe | < 0.15 and |A' 91l .| < 0.09 



2lk\ 



using the OPAL r h v T and tau lifetime [19] measurements and 



A' 31A .| < 0.09 and |A' Q1 J < 0.07 



2\k\ 



using the world average r 



branching ratio and tau lifetime. Previous measurements have quoted 



the 68% (la) confidence limits [22]. Using the OPAL measurements, the 68% (la) confidence limits are 



|A 31fc | < 0.13 and |A 21fe | < 0.06, while the world average results give 68% confidence limits of |A 31fc | < 0.07 
and |A 21 fc| < 0.05. The limit on X' 21k set using this method is competitive with the present best limit 



derived from pion decay to electrons and muons [21]. Limits on A' 31fc have been obtained previously using 
the decay r~ — > ir~iy T [22]. Our new calculation using r~ — > h~v T has several advantages over this method. 

ir~v T ). Additionally, BR(r~ — > h~v T ) is 



7T v T 



Firstly, BR(t —* h u T ) is more precisely known than BR(r 
a directly measured quantity making the limits more experimentally compelling than those derived from 



7T V T . 



The t~ — > h~ 7T°^ T branching ratio is measured to be (25.89 db 0.17 ± 0.29) % and is compared to other 
published branching ratios in figure ||. Some of the branching ratios shown have been corrected by the PDG 
in order to treat the kaon backgrounds in a uniform manner. The t~ — > h~ branching ratio result 

measured in this work agrees well with the PDG average (25.76 ± 0.15%) and the previous measurements. 



The Conserved Vector Current hypothesis [23] can be used to predict the r~ — > ir~ ir°v T branching 
ratio from low energy e + e~ — > tt + tt~ data. A number of predictions for the r _ — ► 7r~ 7r°^ T branching ratio 
have been made JTT| , p4[ . A recent review by Eidelman and Ivanchenko |2l| predicted the branching ratio 
to be (24.25 ± 0.77)%. The r~ — ► h~ -k®v t branching ratio measured in this work should be modified by 
subtracting the r~ — > K" branching ratio O giving a t — > 7r 7r°i/ r branching ratio of (25.73±0.31)%. 

The result measured here is consistent (within two standard deviations) with the CVC prediction. 



The r~ — > h~ > 27r°z^ T branching ratio is measured to be (9.91 ± 0.31 db 0.27) %, in comparison to the 
PDG average for the r~ — > /i~ > 27r°z/ r branching ratio of (10.48 ± 0.35)% p. The PDG number quoted is 
the sum of the average values for the branching ratios of the tau decay to the h~2ir°v T , h~37r v T1 h~4ir u T , 
ir~K u T and K~K u T modes. 



14 



9 Conclusions 



The branching ratios of r — > h v T , r — > h ir°v T and r — > ft. > 27r°i/ r decays have been measured 
with the OPAL detector at LEP. The branching ratios are 

BR(r ^/i^ T ) = (11.98 ±0.13 ±0.16)% 

BR(r _ /T 7r°i/ T ) = (25.89 ±0.17 ±0.29)% 

BR(r- -> hT > 2tt°u t ) = (9.91 ± 0.31 ± 0.27) % , 

where the first error is statistical and the second error is systematic. These new measurements are more 
precise than previous OPAL measurements and supersede those results. 

The r~ — > h~ u T branching ratio measured in this work is found to be in good agreement with previous 
measurements. The ratio of the charged current coupling constants of muons and taus using the t~ — > h~ v T 
branching ratio is found to be g^j g T = (1.018±0.010), consistent with lepton universality. The r~ — ► h~ v T 
branching ratio is used to place limits on supersymmetric .R-parity violating Yukawa couplings. The 
r~ — ► h~ 7r°^ r branching ratio found in this work is in good agreement with the previous results and also 
with the Conserved Vector Current prediction. Finally, the r~ —> h~ > 2i: v T branching ratio measured 
in this work is found to be consistent with the current PDG world average. 
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Figure 1: The distribution of the number of neutral clusters per jet for one-prong selected jets. The 
distribution is shown in both a linear vertical scale (left plot) and a logarithmic vertical scale (right plot). 




Figure 2: The distribution of the energy of clusters identified as radiative photons. 
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Figure 3: (a) The invariant mass distribution of two neutral clusters in data and Monte Carlo for jets 
with two neutral clusters with energy between 0.5 and 9.0 GeV. (b) The invariant mass distribution for 
the neutral cluster pairs selected by the 7r° algorithm. 
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Figure 4: (a) The distribution of the number of 7r°'s in each one-prong jet identified by the 7r° finding 
algorithm, (b) The distribution of the types of 7r°'s identified in one-prong jets. 
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Figure 5: The energy distribution of reconstructed jt s shown in both linear and logarithmic scales, 
(a) and (b) show the energy of 7r°'s in r~ — > h~ 7r°u T jets, (c) and (d) show the energy of 7r°'s in 
r~ — ► h~ > 27t°ij t jets. The structure observed at 9 GeV corresponds to the one-cluster tt° threshold. 
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Figure 6: The invariant mass distribution of jets identified as r~ — > h~ 7r°z/ T (a) and r~ — ► /i _ > 2n°v T 
(b). The mass cuts are indicated on the plots. In figure (a) the diagonally hatched area is t~ — > h~ > 2ir°v T 
jets, the cross hatched area is t~ — > i/ T jets and the dark shaded area is other background. In figure 
(b) the diagonally hatched area is r~ — > /i~ tt°v t jets and the cross hatched area is other background. 




Figure 7: The dE/dx distribution used to measure the r e u c v T background in the r — > /i" 
sample. A background-dominated region is created by selecting jets with 9.25 < d-E/dx < 12.0 keV/cm 
(indicated on the histogram). The cross hatched region corresponds to the r~ — > e~!7 c zy r background 
events. 
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Figure 8: The branching ratios for the r h u T and r — > h~ decay modes for this work are 

compared with previous measurements. The solid band in each plot is the PDG average for that channel. 
The PDG average for the r~ — > h~ v T branching ratio does not include the CLEO measurement while the 
PDG average for the r~ — > h~ branching ratio uses all of the measurements. The error bars shown 

include both systematic and statistical uncertainties. 
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Figure 9: The tau lifetime is plotted as a function of r~ — > h~ u T branching ratio. The shaded band 
shows the Standard Model prediction assuming lepton universality and its width reflects the uncertainty 
associated with the tau mass. The solid point uses the r~ — > h~ u T branching ratio obtained in this 
measurement and the OPAL tau lifetime. The open circle uses the world average for both the t~ — > h~ v T 
branching ratio and tau lifetime. 
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